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We propose a new signature for weakly interacting massive particle (WIMP) dark matter, a spectral 
feature in the diffuse extragalactic gamma-ray radiation. This feature, a sudden drop of the gamma- 
ray intensity at an energy corresponding to the WIMP mass, comes from the asymmetric distortion 
of the line due to WIMP annihilation into two gamma-rays caused by the cosmological redshift. 
Unlike other proposed searches for a line signal, this method is not very sensitive to the exact dark 
matter density distribution in halos and subhalos. 



PACS numbers: 95.35.+d; 14.80.Ly; 98.70.Rz 

It has been known since long that particle dark mat- 
ter annihilations may produce an observable gamma-ray 
line @-§. One of the prime particle dark matter candi- 
dates is a WIMP (Weakly Interacting Massive Particle), 
of which the supersymmetric neutralino is a favourite ex- 
ample. WIMP annihilation into 77 and Z7 would give 
monochromatic gamma rays with an energy equal to (or 
close to) the WIMP mass jlJffi . Since these gamma rays 
are monochromatic and have high energy they would con- 
stitute a spectacular signature of annihilating dark mat- 
ter. 

There has been a rapid development of the understand- 
ing of how structure forms in the Universe. In the current 
model for structure formation, the ACDM model, most 
of the matter is in the form of non-relativistic cold dark 
matter (CDM), but with a contribution to the energy 
density also from a cosmological constant (A) . As shown 
by detailed A^-body simulations (see, e.g., p||8| and ref- 
erences therein), in such a picture large structures form 
by the successive merging of small substructures, with 
smaller objects generally being denser. The A^-body sim- 
ulations also show that the dark matter density profile in 
clusters of galaxies and in single galaxies develops a steep 
cusp near the center, pcDM{r) ~ r _Q with a ranging 
from 1 § to 1.5 @. 

At present, it is not clear whether these A^-body pre- 
dictions are in agreement or not with available data. On 
large scales, this scenario gives excellent agreement with 
observations, see, e.g., the prediction of the Lyman-a ab- 
sorption lines at high redshifts On smaller scales, 
one of the main puzzles is how to properly include bary- 
onic matter. For instance, it appears that the contra- 
diction between the number of satellites found in the N- 
body simulation of a halo with the size of the Milky Way 
and the number of those observed may be explained by 
plausible mechanisms which make most small subhalos 
dark It is less clear how to get agreement between 
the measured rotation curves of dwarf and low surface 
brightness galaxies and those found in ACDM simula- 
tions (see |p| for a recent review). 



Here we will take the view that the ACDM picture 
is basically correct and that structure forms hierarchi- 
cally, with the number density of halos of mass M being 
distributed as dN/dM oc Af"* 3 with /3 ~ 1.9 - 2, as pre- 
dicted by Press-Schechter theory j[4| and also verified in 
A^-body simulations. Furthermore, we will use that the 
concentration of halos grows with decreasing mass. 

Previous analyses (e.g., JIe|-|2"c|1) have focused on the 
dark matter gamma-ray signals from the Galactic center 
and the halo of the Milky Way or isolated nearby galax- 
ies and satellites; in these cases the actual dark matter 
distribution within halos plays a crucial role for the ob- 
servability. The presence of substructures, as well as of 
central cusps, increases the detected rates |IJ|,|l7] |l9| , but 
still it may be difficult to resolve such individual sources 
(see 111] ) . We now show that the integrated signal of un- 
resolved cosmological dark matter halos gives a potential 
detection method which is more robust to changes of the 
details of how the dark matter is distributed locally. 

We consider the lightest neutralino, \, of the MSSM 
(the Minimal Supersymmetric Standard Model) as our 
template particle. The mass range is from around 50 
GeV up to several TeV (see (2l]] for a recent review). 
We start with the (unrealistic) case of all the dark mat- 
ter being smoothly distributed at all redshifts, and then 
modify the results by introducing structure. The comov- 
ing number density n c of neutralinos, after decoupling 
from chemical equilibrium ("freeze-out") at very large 



temperatures (T 



c /20) is depleted slightly due to 



self-annihilations, governed by the Boltzmann equation 
h c = —(<jv)(1 + z) 3 n^,, where (av) is the thermally- 
averaged annihilation rate, which, to an excellent approx- 
imation after freeze-out, is velocity independent, since 
the neutralinos move non-relativistically. 

Each of the \ particles that disappears will give rise to 
Nj photons on the average, with an energy distribution 
in the rest frame of the annihilation 
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where the first term gives the average continuum gamma 
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ray distribution per annihilating \ an d the second term 
is the 77 line contribution, with fr 77 being the branch- 
ing ratio into 77 (in this discussion, we neglect the Z7 
channel |0j). 

Gamma-rays observed today with an energy Eg corre- 
spond to an emitted energy E = (1 + z)Eq. If we now 
track, using the Boltzman equation, the number of neu- 
tralinos that have disappeared from redshift z until now, 
and fold in the energy distribution Eq. ([!]), we can esti- 
mate the diffuse extragalactic gamma ray flux. Let Hq be 
the Hubble parameter. Using the relation between time 
and redshift (see , e.g., p2|) d/dt = —H a (l + z)h (z)d/dz 
with h(z) = ,/%(! + z) 3 + fl K (l + z) 2 + tt Al where 
57 a/, CI a and Qk = 1 — &m — £7 a are the present fractions 
of the critical density given by matter, vacuum energy 
and curvature, the rate is 
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where k — (<rv) /H ■ 

The differential spectrum of the number density n 7 of 
photons generated by annihilations is given by: 



dn 7 
dz 



= AL 



dn c 
dz 



dN 1 (E) dn c 
dE dz 



dE. 
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Here, dn c /dz can be computed directly from (g) to excel- 
lent accuracy, replacing the exact solution n c {z) by the 
present number density of neutralinos no on the right 
hand side. 

Approximating fl x ~ Q,m, we obtain no = p x /m x = 
Pcrit^M / m x' wnere Pcrit = 1-06 • 10~ 5 h 2 GeV/cm 3 and 
h is the Hubble parameter in units of 100 km s _1 Mpc -1 , 
and we find that the gamma-ray flux is given by: 
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8.3 • 10~ 14 cm _2 s- 1 sr- 1 GeV" 
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where we defined r 2 6 = (<7i>)/(10 _26 cm 3 s _1 ) and mioo 
the mass in units of 100 GeV. The term e -z / Zmax accounts 
for absorption of gamma-rays along the line of sight. For 
the energies we are interested in, 1 GeV < Eq < 500 
GeV, it is the starlight and (poorly known) infrared back- 
ground radiation which is the limiting factor. An optical 
depth of order unity is reached for a redshift which can 
be approximated by z max (£ , ) ~ 3.3(£ /10 GeV)" 8 , 
which is a simplified representation of the results of 
p3| , p4| . The exponential form is a good approximation 
for small values of z max as in most of our cases. The 
upper limit of integration is z up = m x / Eq — 1, since the 
maximum rest frame energy of a photon in an annihi- 
lation event is E = m x . The gamma line contribution 
to (Q) is particularly simple, just picking out the inte- 
grand at z + 1 = m x /i?o; it has the very distinctive and 



potentially observable signature of being asymmetrically 
smeared to lower energies (due to the redshift) and of 
suddenly dropping just above m x . The continuum emis- 
sion will produce a less conspicuous feature, a smooth 
"bump" below one tenth of the neutralino mass which 
may be more difficult to detect. To give an example of 
the results, we performed a large scan of the MSSM pa- 
rameter space obtained with the DarkSUSY package [^5| . 
Models with large 77 rates ((<7u)2 7 ^ 10~ 29 cm 3 s _1 ) ex- 
ist in all the mass range from m x — 70 GeV to several 
TeV. Consider a high-rate model with m x = 86 GeV, 
T 2 6 ~ 6, 6 77 ~ 3 • 10 -3 , in the "concordance" cosmology 
fl M = 0.3, Q A = 0.7, h = 0.65 (see Ref. @ for the 
full range of predicted fluxes in the MSSM). The con- 
tinuous gamma-ray rest frame energy distribution per 
annihilating particle comes mainly from hadronization 
and decay of tt°s and is conveniently parametrized as 
dN cont (E)/dE = {0A2/m x )e- 8x /{x 1 - 5 + 0.00014) where 
x = E/m x . The resulting flux near 86 GeV is around 
10 -15 cm~ 2 s~ 1 sr _1 GeV~ 1 , some five orders of magni- 
tude below the diffuse extragalactic flux extrapolated 
from the Energetic Gamma-Ray Experiment Telescope 
(EGRET) measurements |fj. 

We now include the important effects of structure for- 
mation. Consider a halo of mass M whose radial density 
profile can be generically described by |2^] pdm{t) — 
p' DM f (r/a), with p' DM being a characteristic density and 
a a length scale. These are found in V-body simula- 
tions not to be independent parameters, as smaller halos 
are associated with higher densities. We assume that 
the halo of mass M accreted from a spherical volume 
of radius Rm, determined by requiring that the average 
cosmological density times that volume is equal to M, 
po ■ <hrR 3 M /3 = M (with p Q - 1.3 ■ 10~ 6 GeV/cm 3 ). The 
increase of average squared overdensity per halo (which 
is what enters the annihilation rate) is given by: 
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where /„ = J B ' M ^ a y 2 dy(f(y)) n . Here the dependence on 
the upper limit of integration is rather weak, while for 
the lower limit of I2, in very singular profiles, like the 
Moore profile, a cut-off has to be imposed due to rapid 
self-annihilations near the center |20[ . 

Computing I2/I1 numerically, and using values of 
p'riM / Po as determined for Milky Way size halos from 
|g , we find values of A 2 of 2.3 • 10 5 for the Moore profile, 
1.5 • 10 4 for the Navarro-Frenk- White (NFW) profile [|, 
and 7- 10 3 for a cored, modified isothermal profile (modi- 
fied such that the density falls as 1/r 3 at large radii p0] ) . 
The flux ratios, 30 : 2 : 1 for these three models should be 
compared with the ratios 1000 : 100 : 1 obtained within 
a 5-degree cone encompassing the galactic center ]2(| . 

We also take into account that the number density of 
halos is scaling like ~ 1/M 2 , and that small-mass halos 
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are denser. Again, we resort to the highest-resolution 
iV-body simulations available to date. Fitting the con- 
centration parameter of Moore-type halos by p(J c ~ 
100 (M vir //i _1 M Q )~ 084 , we find to a good approxima- 
tion A 2 - 2 • 10 5 M 12 ' 22 , where M i2 is the halo mass 
in units of 10 12 solar masses. This means that the to- 
tal flux from a halo of mass M scales as M 78 . Since 
the number density of halos goes as M~ 2 , the fraction 
of flux coming from halos of mass M scales as M -122 . 
Thus the gamma-ray flux will dominantly come from the 
smallest CDM halos. In simulations, substructure has 
been found on all scales (being limited only by numerical 
resolution). For very small dark matter clumps, however, 
no gain in overdensity is expected, since once the matter 
power spectrum enters the fc~ 4 region a constant density 
is predicted |jg|| . We conservatively set 10 5 — 10 6 Mq as 
the minimal scale. 

Finally, regarding redshift dependences, we have as- 
sumed a constant enhancement factor A 2 out to z ~ 1, 
and have arbitrarily imposed quadratic growth in the en- 
hancement factor from z ~ 10 to the fully developed 
epoch z = 1. (The observable feature is not sensitive to 
this assumption.) Furthermore, in Eq. (^) we make the 
replacement (1 + z) 3 — > 1 [^7|, reflecting the fact that 
the we are now dealing with a clumped, rather than a 
smooth distribution with density scaling ~ (1 + z) . 

We thus arrive at the following expression for the flux 
including structure formation 
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We find that the flux from small Moore-like halo struc- 
ture is enhanced by roughly a factor 2 • 10 6 compared to 
the smooth case, giving observability for the same anni- 
hilation parameters as used above. In Fig. [I], we show 
the results for the same 86 GeV MSSM model as in the 
example above, and a model of 166 GeV mass, = 59, 
& 77 = 1.2 ■ 10~ 4 . 

Several remarks are in order here. 

(a) The flux calculated here is made up by the inte- 
grated effect of a large number of halos. One could also 
detect nearby clumps in the Galactic halo [^,[[9). This 
would in particular make the line even more visible. 

(b) Contrary to previous expectations (e.g., f2£|), we 
do not find it excluded that a large part of the measured 
extragalactic gamma-ray flux around and above 10 GeV 
may originate from WIMP annihilations. Our more op- 
timistic conclusion rely primarily on three factors: (1) 
The improved understanding of CDM structure forma- 
tion, through semianalytic and numerical methods show- 
ing that less massive halos are denser. (2) The possibility 
that the density profile of each subhalo can be steep, in- 
creasing further the emissivity. (3) The observability is 
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FIG. 1. The predicted diffuse 7-ray flux, from cosmic anni- 
hilations into continuum gamma-rays, and a gamma-ray line. 
The redshifted line gives the conspicuous feature at the high- 
est energies. Shown are cosmic annihilation of 86 GeV (solid 
line) and 166 GeV (dotted line) neutralinos. A Moore den- 
sity profile for the halo substructure has been assumed. The 
EGRET data |2^] on the extragalactic flux are the data points 
with error bars shown. 



helped by the existence of a sharp gamma-ray line in the 
annihilation spectrum and the expected opacity of the 
universe to this signal at the energies considered here, 
which gives a distinct feature with no known astrophys- 
ical background. None of these three ingredients were 
present in the older analyses. 

(c) An advantage of looking for a spectral feature in 
the diffuse extragalactic gamma-ray background is that 
regions of the sky can be used where foreground con- 
tamination is as low as possible. Previous search strate- 
gies, focusing on the Milky Way center or the center of 
other nearby galaxies, are hampered by the relatively 
high gamma-ray flux from other sources along these lines 
of sight. 

(d) We note that the GLAST satellite [§§ will cover 
the energy range up to 300 GeV with unprecedented pre- 
cision. It is possible that most of the measured EGRET 
flux comes from Active Galactic Nuclei (AGN). Then, 
the absorption on the infrared background will be se- 
vere above 80 GeV for all but the very nearest AGNs. 
The true diffuse background above 80 GeV may there- 
fore be one or two orders of magnitude smaller than 
the EGRET extrapolation ~ E~ 21 , with corresponding 
higher chances of detecting a line signature. Also note 
that the flux from the 27 line and from continuum gam- 
mas in neutralino annihilation are not strongly correlated 
(except in the case of very pure higgsinos |16}| ), so it may 
well be that the line is visible but not the "bump" , and 
vice versa. 

(e) We notice that the spectral features of Fig. |l| would 
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appear with high significance, for the line it is of the order 
of 10cx for a 5-year exposure with GLAST assuming the 
background extrapolation ~ E~ 21 given in pfJ], which 
becomes 20a if, in analogy with the spectral slope of other 
cosmic rays, the background above 10 GeV drops like 
E~ 2J instead. Even if halos have less singular density 
profiles, such as the cored isothermal sphere, a detectable 
signature may appear. Alternatively, for the Moore and 
NFW profiles, models with smaller gamma-ray branching 
ratios and/or higher masses may be probed. 

(f) The relatively broad spectral feature caused by the 
redshifted gamma-ray line alleviates the demand for very 
high energy resolution of the detecting instrument put by 
the previous suggestions of searching for the line in the 
local Galactic neighbourhood. We find that an energy 
resolution of 10-20 % is in fact adequate. For GLAST, 
this means that the effective area can in fact be larger 
than we have assumed (since a smaller requirement on 
the energy resolution means that events from a larger 
part of the field of view of the detector can be used). 

(g) It has to be reminded that the strength of the 
gamma-ray line can be much lower than the examples we 
chose for Fig. [|, in which case a detection would be cor- 
respondingly more difficult. However, in many of those 
cases the continuum signal may be large. In particular, 
we find that models compatible with the recent results 
on the muon g — 2 Q according to the suggested contri- 
butions of "new physics" (such as supersymmetry) of pl| ] 
all give high values for the continuum gamma-ray flux. 

(h) Although the clumpy structure is in the dark mat- 
ter distribution, and may not have optical counterparts, 
the rate of the annihilation should scale with the overall 
matter distribution in the nearby Universe. Thus, once 
mass maps from, e.g., the Sloan Digital Sky Survey are 
available, a cross-correlation analysis with the gamma- 
ray maps should enable a further reduction of Galactic 
gamma-ray foregrounds. 

To conclude, we have suggested a new possible signa- 
ture for dark matter detection, which employs the clus- 
tering properties of Cold Dark Matter as they emerge 
from A-body simulations. The upcoming GLAST satel- 
lite, and possibly ground based detectors, may have a 
good chance at detecting the characteristic features, in 
particular the rise and sudden drop in gamma-ray flux 
around the WIMP mass produced by annihilation into 
two photons. 
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